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• • ■ Charged kaon production has been measured in Si+Al and Si+Au collisions at 14. 6A GeV/c, and 

Au+Au collisions at 11.1 A GeV/c by Experiments 859 and 866 (the E-802 Collaboration) at the 
Jv^ , BNL AGS. Invariant transverse mass spectra and rapidity distributions for both K + and K~ are 

presented. The centrality dependence of rapidity-integrated kaon yields is studied. Strangeness 
enhancement is observed as an increase in the slope of the kaon yield with the total number of 
participants as well as the yield per participant. The enhancement starts with peripheral Si+Al and 
Si+Au collisions (relative to N+N) and appears to saturate for a moderate number of participating 
nucleons in Si+Au collisions. It is also observed to increase slowly with centrality in Au+Au 
collisions, to a level in the most central Au+Au collisions that is greater than that found in central 
Si+A collisions. The enhancement factor for K + production are 3.0 ± 0.2(stat.) ± 0.4(syst.) and 
4.0 ± 0.3(stat.) ± 0.5(syst.), respectively, for the most central 7% Si+Au collisions and the most 
central 4% Au+Au collisions relative to N+N at the correponding beam energy. 

PACS number(s): 25.75.-q 
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I. INTRODUCTION 

Strange particle production has long been considered 
as a useful experimental probe of heavy-ion collisions be- 
cause strange particle yields are thought to be sensitive to 
physics resulting from the high baryon densities attained 



in heavy-ion collisions at AGS energies In particular, 
quark-gluon plasma formation is expected to result in en- 
hanced strangeness production [^|-|6). Such an enhance- 
ment was observed by BNL Experiment 802 in Si-induced 
collisions as an increase in the K + yield in central colli- 
sions relative to the pion yield [0,^). However, secondary 
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interactions among hadrons have been proposed to ex- 
plain the observed strangeness enhancement, and cas- 
cade models implementing these secondary interactions 
can approximately reproduce the kaon yields measured 
in Si+A collisions [p|-|l2|]. Because there are, in principle, 
many mechanisms that can increase the yield of strange 
particles in heavy-ion collisions, the process(es) respon- 
sible for the observed strangeness enhancement have not 
experimentally been clearly identified. However, with the 
present high statistics Si+A and Au+Au data, the sys- 
tematics of strangeness production can now be studied 
as a function of the system size from very small systems 
to very large systems, thus, providing more stringent 
constraints on theoretical explanations of strangeness en- 
hancement. 

In this paper, semi-inclusive K + and K~ spectra mea- 
sured by E-859 in Si+Al and Si+Au collisions at 14.6^4 
GeV/c and by E-866 in Au+Au collisions at 11. 1A 
GeV/c at the BNL AGS are presented. Both the Si+A 
and Au+Au measurements were made by the E-802 col- 
laboration using the Henry Higgins spectrometer Jl3| . 
Previous measurements of kaon production in Si+A col- 
lisions have been presented by E-802 • Kaon data 
from Au+Au collisions at slightly higher beam energy 
are presented by E-866 in ATE). 

This paper is organized in the following way. First, 
descriptions of the experimental apparatus and data col- 
lection are given in section II. Then, the data analysis 
is described in section III. Results on kaon spectra are 
described in section IV, followed by discussions on the 
centrality dependence of the kaon production rates in 
section V. Finally, conclusions are drawn in section VI. 



II. EXPERIMENTAL APPARATUS AND DATA 
COLLECTION 

E-859 was an upgraded experiment from E-802 at the 
BNL AGS. The upgrade consisted of an addition of two 
multi-wire proportional chambers for triggering purpose 
to the E-802 spectrometer between the Henry Hig- 
gins dipole magnet and the time-of-fiight wall, and an 
implementation of a level-2 trigger 01J . The level-2 trig- 
ger linked hits from one of the trigger chambers and the 
time-of-fiight wall by straight lines to be verified on the 
other trigger chamber. Verified combinations were as- 
sumed to originate at the target and particle momenta 
were obtained. Particle masses were then obtained from 
the momenta and the time-of- flight, and were used for 
event selection. Event selection was achieved on-line by 
looking up a series of tables pre-loaded in LeCroy CA- 
MAC modules. For most of the data presented here, 
the level-2 trigger rejected events not containing a K + , 
K~ , or p. K + and K~ (or p) were defined as particles 
for which the calculated masses fall in the range of 0.3 
< m < 0.7 GeV/c 2 and 0.3 < to < 1.5 GeV/c 2 , re- 
spectively. Events were also taken with an override on 



the level-2 trigger decision to check biases in the level-2 
trigger selection criteria. 

E-866 was an upgraded experiment from E-859 de- 
signed for measurements of the Au+Au collisions. The 
main upgrade was the addition of a spectrometer for for- 
ward angle measurements (TjJTTjJl^]. The upgrade rele- 
vant to this paper included rebuilding the global detec- 
tors for event characterization, replacing the drift cham- 
ber closest to the target, and an addition of two highly 
segmented multi-wire chambers [jlj],[l^[l9| in front of the 
Henry Higgins magnet. The addition of the wire cham- 
bers was driven by the expected high multiplicity of 
central Au+Au events and the need to eliminate back- 
ground tracks. Both wire chambers had drift-time read- 
out. Typical resolution of 200 fim was achieved for both 
chambers. 

The E-859 data were taken during the 1991 and 1992 
Si-beam runs at the AGS, and were an improvement over 
the original E-802 data mainly due to increased statistics 
afforded by the level-2 trigger and the improved accep- 
tance at low momentum. The data were obtained using 
a beam of nominal intensity 10 6 ions/s with Al and Au 
targets of thickness 817 and 944 mg/cm 2 , corresponding 
to a 3% and an 1% Si interaction rate, respectively. The 
data were obtained using 10° steps in spectrometer angle 
covering the range 5-58° and using fields of ±2, ±4 KG 
in the spectrometer mag net [|0U] 



21 



The E-866 data were taken during the Fall 1993 Au- 
beam run at the AGS. The typical beam intensity was 
10 5 ions/s. For the data presented in this paper, a 
975 mg/cm 2 thick Au target was used, corresponding 
to a 1.5% Au interaction rate. The data were obtained 
using 5° steps in spectrometer angle covering the range 
21-58° and using fields of ±2 KG in the spectrometer 
magnet. 

Both data sets were obtained using a global level-1 
trigger that selected interactions by vetoing on the pres- 
ence of a beam-charge ion in the Bull's-eye trigger |t3 22 
downstream of the target, with the additional require- 
ment in the spectrometer of at least one hit on the time- 
of-fiight wall and a hit on either of the two trigger cham- 
bers. A plastic scintillator was used in the Bull's-eye for 
the Si+A data taking with a threshold, corresponding to 
a cut on the remaining charge in the projectile of ap- 
proximately Z < 12.6. For the Au+Au data taking, the 
Bull's-eye was re-designed and rebuilt as a Cerenkov ra- 
diator. The new Bull's-eye had roughly the same geomet- 
rical area as the old one, and used quartz as its radiator. 
A Bull's-eye threshold, corresponding to a cut on the re- 
maining charge in the projectile of approximately Z < 73, 
was used. The cross-section for events satisfying this cut 
was measured to be 5.35 barn. (The total cross-section of 
Au+Au reaction is 6.8 barn |ljj2^].) Further filtering of 
data was performed by the level-2 trigger, when enabled, 
to enhance the sample of particular particle species in 
the data as discussed above. 
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III. DATA ANALYSIS 

In the offline analysis, cuts were applied to the beam 
counter and Bull's-eye data to remove pile-up and to re- 
fine the online Bull's-eye cut. The refined Bull's-eye cut 
corresponded to a remaining charge of Z < 1 1 in the pro- 
jectile for the Si+A data, and Z < 69 in the projectile 
for the Au+Au data, respectively. The cross-section for 
events passing this cut was 1.37 barn for Si+Al, 3.68 barn 
for Si+Au, and 5.23 barn for Au+Au. For all events sat- 
isfying the pile-up and Bull's-eye cuts, tracks were re- 
constructed and identified in the spectrometer. For the 
Si+A data, track segments were reconstructed from hits 
in the drift chambers behind the magnet, and projected 
through the magnetic field to the front of the spectrom- 
eter. The projections were verified by hits in the front 
drift chambers which were refit to give final track pa- 
rameters. For the Au+Au data, the track reconstruction 
behind the magnet was done as for Si+A, and an in- 
dependent track reconstruction was performed in front 
of the magnet jl9| . The front track reconstruction used 
the highly segmented wire chambers to define narrow re- 
gions (on order of a few wire pitches, or 1 cm) for track 
candidates. Front drift chamber hits found in these re- 
gions were then used, together with the wire chamber 
hits with drift time information included, to obtain track 
segments in the front. The track segments in the front 
and in the back were then matched through the mag- 
netic field. A large fraction of background tracks did 
not have both the front and the back segments, therefore 
the matching algorithm effectively rejected background 
tracks. Good tracks were required to project within 2 
cm around the nominal target position, further reducing 
background tracks. 

The Au+Au particle identification was performed us- 
ing data from the time-of-flight wall, since the spectrom- 
eter was operated at larger angles where particle mo- 
menta are typically small. A particle was identified as 
a kaon if its measured inverse velocity, 1//3, was within 
±3(7 of the expected value and over 3cr away from the 
expected values of other particle types. The expected 
values for various particle types were calculated from the 
time of flight, track path length up to the time-of-flight 
wall, and the corresponding particle masses. The typ- 
ical track path length up to the time-of-flight wall was 
660 cm. The a was a combination of time-of-flight reso- 
lution and deviation of track trajectory from straight line 
due to multiple scattering. The time-of-flight resolution 
was measured to be 130 ps. Hence, kaons were identified 
up to a momentum of 1.74 GeV/c. For the Si+A data 
both the time-of-flight wall and the segmented Cerenkov 
counters (l^] were used for particle identification. The 
time-of-flight resolution was measured to be 120 ps for 
these data. Kaons were identified up to a momentum 
of 3 GeV/c. For both data sets, the contamination of 
the K + and K~ samples was found to be below 1% and 
3%, respectively, averaged over the spectrometer accep- 



tance JTs| ] . 

In order to study the centrality dependence of kaon 
production, the data were divided into bins in centrality. 
For the Si+A data the centrality selection was made us- 
ing the E-859 Target Multiplicity Array (TMA), and for 
the Au+Au data the selection was made using the zero- 
degree calorimeter (ZCAL). The Si+A cuts were made in 
multiplicity because of a rate-dependent worsening of the 
zero-degree calorimeter resolution during the high beam 
intensity E-859 run pi]] . The zero-degree calorimeter 
data was used to determine the correlation between mul- 
tiplicity and zero-degree energy, E, calculated for events 
in each multiplicity bin. The results are shown in Table |. 
The centrality cuts applied to the data-sets are also listed 
in Table [| along with the cross-sections for interactions 
satisfying the cuts. 

Due to the large amount of radiation damage from 
the high beam intensity run in E-859, the zero-degree 
calorimeter was rebuilt for the Au+Au run. In the 
Au+Au analysis, the zero-degree energy was calibrated 
run-by-run to account for a gradual degradation in the 
calorimeter response due to radiation damage p9|. The 
total degradation over the whole run period was 5%. The 
zero-degree energy resolution was measured to be 1.48 
GeV 1 / 2 • \f~E. The calibrated and target-out subtracted 
zero-degree energy spectrum for Au+Au events is shown 
in Fig. m. The target-out zero-degree energy was shifted 
slightly lower to compensate for the average energy loss of 
the beam in the target prior to subtraction. The target- 
out interaction rate was 1.1% for Au+Au and 1.0% for 
Si+A. Also shown in Fig. |l| as the dotted lines are the 
centrality cuts used for the Au+Au data. The mean zero- 
degree energy for each bin are determined from the dis- 
tribution and results are shown in Table |. 

Total number of identified K + 's and A~'s, after all 
cuts, were 120K and 44K for Si+Al, 151K and 69K for 
Si+Au, and 28K and 11K for Au+Au. The statistics were 
uniformly distributed in available centrality bins. To ob- 
tain semi-inclusive spectra, A^ + 's and isf - 's from events 
satisfying each centrality cut were binned in transverse 
momentum p± , and rapidity y = | In , where (3 Z is 
the longitudinal velocity of the particle. Invariant differ- 
ential yields, azimuthally averaged, per collision within 
a given centrality bin, were obtained according to the 
formula, 

d 3 N d 2 N d 2 N 

dp 3 27rm±dm±dy 2irp±dp±dy 
1 AN(y, P± ) 1 
A cv t p±Ap ± Ay A(y,p±)' 

where A cvt is the number of events that fell within the 
given centrality bin; the transverse mass, m±, is defined 
as mj_ = \J v\ + m o ; Ap±, Ay are bin sizes in p±, y, re- 
spectively; AN(y,p±) is the number of kaons in the bin; 
and A(y,p±) is the correction factor. For the Au+Au 
data, Ap ± = 50 MeV/c 2 for both K+ and K~; Ay = 0.2 
for K~\ Ay = 0.1 for K + in all bins except the most 
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peripheral where Ay — 0.2. The binning is compara- 
ble to that used in the Si+A data set. Only particles 
with momentum 0.3 < p < 1.74 GeV/c were used in 
the spectra. The correction factor, A(y,p±), accounted 
for the azimuthal acceptance of the spectrometer, parti- 
cle decays, multiple-scattering and hadronic interaction 
losses, track reconstruction and particle identification in- 
efficiency [jigj] . The spectrometer acceptance was calcu- 
lated from the geometry of the detectors. Other effects 
were determined by Monte Carlo simulations. The par- 
ticle decay correction for kaons was found to be roughly 
the same as the theoretical momentum-dependent cor- 
rection because very few kaons that decayed before the 
time-of-flight wall were actually reconstructed. Track re- 
construction efficiency depended on hit multiplicity in 
the chambers, and was 98% with the spectrometer at 44°, 
dropping to 95% at 21° for the Au+Au data ||]. For the 
Si+Au data, it was above 90% with the spectrometer at 
5° and was better at larger spectrometer angle settings. 
For the Si+Al data, it was better than that for the Si+Au 
data at the same spectrometer angle settings. Track loss 
due to hadronic interaction for K + was 1.1%, indepen- 
dent of momentum; for K~ it was 2.2% at high momen- 
tum (> 0.6 GeV/c) and increased to 10% at momentum 
0.3 GeV/c. Track losses due to multiple scattering and 
particle identification inefficiency were combined. The 
net effect was the same for K + and K~, and was 4.5% 
for high momentum kaons (> 1.4 GeV/c), and increased 
to 8.3% at momentum 0.6 GeV/c, and 38% at momen- 
tum 0.3 GeV/c. The level-2 trigger inefficiency was mea- 
sured to be less than 1%, and no correction was applied 
to the spectra. Chamber efficiencies were above 95% for 
all sense wire planes in the chambers. Since not all sense 
wire planes were required for reconstructed tracks, the 
effect of hardware inefficiencies was negligible |L9f . The 
overall systematic error is estimated at 5%-10%~^ 

For the Au+Au data, the spectra are corrected for the 
energy loss of particles in the target and in air. The en- 
ergy loss in the target was calculated on average using 
half thickness of the target and the particle emitting an- 
gle. The energy loss of minimum ionizing particles was 
about 2 MeV. The energy loss of a 0.3 GeV/c momentum 
kaon was 7.5 MeV, corresponding to a momentum shift 
of 5%. 



IV. RESULTS 

Sample sets of kaon transverse mass spectra from 
Si+Au and Au+Au collisions, each for the most central 
bin, are shown in Fig. || and Fig. ||, respectively. All 
semi-inclusive invariant spectra including those shown in 
the figures are well described by an exponential fall-off 
with mi. These spectra are fit to the form 



E- 



d 3 N dN 



dp 3 



-(m ± -m )/T 



dy 2vr(m T + T 2 ) 



(2) 



which implicitly integrates the spectra over the entire m± 
range, and allows the dN/dy values and the correspond- 
ing errors to be directly extracted from the fit. Results 
of these fits are superimposed on the spectra in Fig. || 
and Fig. || as dotted lines. The obtained dN/dy values 
for Si+Au collisions for all four bins in multiplicity mea- 
sured in the TMA are plotted in Fig. |. For Au+Au col- 
lisions, for clarity, only four of the available seven bins in 
zero-degree energy are shown. The error bars include sta- 
tistical errors from the fits and 5% bin-to-bin systematic 
errors. The systematic errors were estimated by com- 
paring yields in the overlap region between spectrometer 
settings and by studying the sensitivity of the spectra to 
the m± intervals used in the fits. In addition, there was 
10% overall systematic error in normalization. 

For the purposes of this analysis, rapidity integrated 
yields were chosen to characterize the kaon production 
rates. The Si+A dN/dy spectra, were fit to Gaussians 
and integrated to get the total kaon yields. The re- 
stricted rapidity coverage of the Au+Au data can be 
compensated for by exploiting the required symmetry 
of the Au+Au system about the nucleon-nucleon cen- 
ter of mass rapidity (yNN = 1-59). If we assume that the 
Au+Au kaon spectra are well-behaved at mid-rapidity 
(as they are in [15]), then the data are sufficient to ob- 
tain an estimate of the total kaon yield for collisions in 
each centrality bin. 

The Gaussian fits which account for the unmeasured 
tails of the rapidity distributions for both data sets are 
shown in Fig. ^ as dotted curves, and were integrated 
over all rapidity to estimate the total kaon yields, the 
results of which are summarized in Table Q. The quoted 
errors on the kaon yields and dN/dy Gaussian cr's are 
statistical only. Systematic errors on the Si+A yields 
were estimated to be 10%. Comparisons with data from 
the E-866 forward spectrometer [ ^5|Jl^ ] indicate that the 
assumption regarding the shape of the Au+Au kaon spec- 
tra at mid-rapidity is reasonable, leading to an estimate 
that the extrapolation to mid-rapidity introduces a few 
percent uncertainty into the present Au+Au total yields. 
The fractions of the estimated total K + and K~ yields 
in the extrapolated tails are typically 30% and 20%, re- 
spectively. The systematic error on the extracted total 
kaon yields in Au+Au is estimated to be of order 10% 
by comparing the extracted yields to alternative fits that 
have different shapes for the tails of the distributions. 
Using the sum of the measured dN/ dy instead of the es- 
timated total yields for Au+Au collisions does not alter 
the dependence of the yield on collision centrality except 
for a scale factor. In order from most peripheral to most 
central collision, the sum of the K + dN/dy measured in 
rapidity range 0.5-1.3 are 0.29, 1.08, 2.37, 3.56, 4.54, 5.38 
and 6.05; the sum of the K~ dN/dy measured in rapidity 
range 0.4-1.4 are 0.076, 0.23, 0.52, 0.74, 0.94, 1.09 and 
1.28. 

For Au+Au collisions, the integrated yields from 
the present data have been compared with those de- 
scribed in Jl5|] taken at a slightly higher beam energy 
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(11.6 A-GeV/c). The K+ and K~ yields of the two 
measurements were compared at the 5 centrality cuts 
from p~5| . The ratios of the yield at the higher en- 
ergy to that found in the present work are 1.12 ± 0.03 
and 1.19 ± 0.04 for K + and K~, respectively, where the 
quoted errors are statistical only. From a parameterized 
fit to p+p data |p5| , the expected ratios are 1.06 and 
1.12 for K + and K~ , respectively. The differences be- 
tween the two data sets, when corrected for the difference 
in beam energy, are well within the estimated overall sys- 
tematic error of 10% in the measurements. The results 
from the 11.1^4 GeV/c beam energy were used in the 
present study because this set contained higher quality 
data for kaon production in peripheral collisions. 

V. DISCUSSION 

In order to discuss the centrality dependence of the 
kaon yields, it would be desirable to plot these yields 
against a more intuitive variable than the zero-degree 
energy or multiplicity. In the Wounded Nucleon Model, 
the particle yields in heavy-ion collisions are expected to 
increase roughly in proportion to the number of projectile 
and/or total participants (participating nucleons) p6| , p7[ . 
Results from previous Si+A measurements || suggest 
that kaon yields within a limited fixed rapidity window 
increase faster than linearly with the number of projec- 
tile participants for Si+Cu and Si+Au, but not for Si+Al 
collisions. However, this procedure treats the target and 
projectile nucleons very differently for asymmetric colli- 
sion systems. Consequently, for this work we have chosen 
to use N p , the total number of participants to character- 
ize collision centrality. 

The energy measured in the zero-degree calorimeter 
results predominantly from projectile spectators, so the 
average number of projectile participants for a given av- 
erage zero-degree energy can be estimated using the for- 
mula 

^pp =A p x(l- -J!—) , (3) 

\ I-beam J 

where E is taken as the average zero-degree energy in the 
centrality bin, A p is the number of nucleons in the projec- 
tile (28 for Si and 197 for Au), and Ef, eam is the total ki- 
netic energy of the beam, 383.4 GeV for Si beam and 2000 
GeV for Au beam ^8|, respectively. For the symmetric 
Au+Au and nearly symmetric Si+Al collisions, the total 
number of participants is, on average, simply twice the 
number of projectile participants. For Si+Au collisions, 
the total number of participants has been estimated by 
using nuclear geometry with Wood-Saxon density distri- 
bution and the Glauber model to relate the average num- 
ber of projectile participants, N pp , to the average number 
of total participants, N p , in each centrality bin. Results 
of this calculation are consistent with those obtained in 
previous analyses [ p~4|j29| ] . The estimated total number of 



participants for the centrality bins of the three systems 
are listed in Table H. The statistical errors on N p are 
negligible due to the large number of minimum bias in- 
teraction triggers that were accumulated for each of the 
three systems. The systematic errors on N p for the Si+A 
data were estimated to be 10% by using two independent 
analyses P,pl|. For the Au+Au data the systematic er- 
rors on N p were estimated to be 8% |is|| . The target-out 
interaction rate was subtracted from the ZCAL distribu- 
tion before A^p was calculated adding negligibly to the 
estimate of the overall systematic error. The correlation 
between the ZCAL energy and another measure of cen- 
trality, the transverse energy, has been analyzed by the 
E802 Collaboration and is discussed in Refs. p9|-|3"l|]. 

The total K + and K~ yields from Si+Al, Si+Au and 
Au+Au collisions are plotted in Fig. | as a function of 
Np. One immediate striking feature of this plot is that 
the measured kaon yields from Si+A and Au+Au colli- 
sions have very different dependence on N p . The kaon 
yields from Si+A increase approximately linearly with 
Np with data from the two targets following approxi- 
mately the same line. The kaon yields from Au+Au col- 
lisions, in contrast to those from Si+A, increase faster 
than linearly with N p as can be clearly seen in Fig. 0. 
On the other hand, in each of the three systems, the 
and K ~ yields have the same dependence on N p , differ- 
ing by only a constant scale factor. This is illustrated in 
Fig. H where the ratios of K + to K~ yield are plotted 
against N p . The ratios show a rather weak dependence 
on N p , and are systematically different for different sys- 
tems. (Note that Si+Al and Si+Au are at the same beam 
energy, and Au+Au is at a lower beam energy.) The sim- 
ilarity between K + and K~ is surprising given that their 
production is expected to proceed through very different 
mechanisms, resulting in different dN/dy shapes. For 
all three systems discussed, the K~ dN/dy distributions 
are systematically narrower than the K +, s, as seen from 
Table | 

It is of interest to examine the power-law behavior of 
the kaon yield versus N p : Nk oc N p . Fits to the Au+Au 
data yield the results a = 1.54 ± 0.05 and 1.52 ± 0.07, 
respectively, for K + and K~ . The power factors for K + 
and K~ are essentially same, consistent with the constant 
K + to K~ ratio, and are considerably larger than 1. 

Although the Au+Au kaon yields increase faster than 
linearly with N p , the absolute yields fall significantly be- 
low the Si+A yields for the same number of participants 
(N p ~ 50), which is further shown in Fig. [7| where the 
kaon yield per participant is plotted against A^ p . In order 
to compare these two data sets, the difference in beam 
energy for the two projectiles must be taken into account. 
(The nucleon pair center of mass energy y/s is 5.39 GeV 
for Si+A, and 4.74 GeV for Au+Au). The energy de- 
pendence of kaon production in heavy ion collisions is 
unknown; data from elementary collisions are used to es- 
timate the effect. K + and JV yields in p+p collision at 
12 GeV/c were measured || to be 0.048 and 0.0075, re- 
spectively. Based on a parameterization of p+p data at 
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various energies (25), K + and K yields in p+p should 
be scaled, respectively, by 1.27 and 1.53 from beam mo- 
mentum 12 GeV/c to 14.6 GeV/c, and by 0.90 and 0.82 
from beam momentum 12 GeV/c to 11.1 GeV/c. There- 
fore, K + and K~ yields in p+p at 14.6 GeV/c can be 
estimated as 0.061 and 0.011, and at 11.1 GeV/c can be 
estimated as 0.043 and 0.0062. (Typical errors on these 
yields are 15%.) Since there are significant differences in 
kaon production in p+p, p+n and n+n interactions p3[ , 
the iso-spin averaged kaon yields should be used. Us- 
ing the method described in fl33[ and data measured in 
p+p at 12 GeV/c f^jMj , the ratios of cross sections for 
producing K + in these three elementary interactions can 
be estimated as p+p : p+n : n+n = 1 : 0.81 : 0.62, 
and those for producing K~ can be estimated as p+p 
: p+n : n+n = 1 : 1.19 : 1.42. (See appendix of JT|] 
for details.) Assuming that these ratios do not vary with 
beam momentum between 11.1 GeV/c and 14.6 GeV/c, 
and assuming that the heavy-ion collisions are simple su- 
perpositions of nucleon-nucleon (N+N) collisions at the 
full beam energy, the K + and K~ yields per N+N in- 
teraction can be estimated, respectively, as 0.049 and 
0.013 for Si+A f|| at 14.6 GeV/c and 0.033 and 0.0077 
for Au+Au at 11.1 GeV/c. The typical errors of these 
estimates are on the order of 15%. These results are 
shown in Fig. at N p = 2 as the filled square for the 
Si beam energy and the filled circle for the Au beam 
energy. Based on these results, the Si+A K + and K~ 
yields should be approximately 40% and 90% larger than 
the Au+Au data if the heavy-ion collisions were simple 
superpositions of N+N collisions at appropriate beam en- 
ergy. However, even after compensating for the energy 
difference the Si+A yields are still systematically 100% 
larger than the Au+Au yields for the same N p . 

The above observed differences among the three sys- 
tems could come from the different initial nuclear geom- 
etry. At the same N p , the initial geometry of the three 
systems is by no means the same. For the same N p (^ 50) 
at which the Au+Au collision is peripheral, whereas the 
Si+A collisions are not, the overlap region in Au+Au col- 
lision is less dense than that in Si+Al or Si+Au collision 
due to the dilute surfaces of the nuclei. Therefore, one 
would naively expect lower kaon yield in Au+Au than in 
Si+A collisions at same beam energy, due to fewer suc- 
cessive N+N collisions per participant in Au+Au than 
in Si+A collisions. However, detailed study of the aver- 
age number of successive N+N collisions per participant, 
based on the Glauber model and the inelastic N+N cross- 
section 30 mb, does not give quantitative support to this 
supposition. 

For both Si+Al and Au+Au collisions of all centrali- 
ties, and for peripheral Si+Au collisions, the kaon yields 
per participant increase with N p , as seen from Fig. |t]. 
In non-peripheral Si+Au collisions, however, the kaon 
yield per participant saturates at the value observed in 
the most central Si+Al data. The saturation cannot be 
explained by the simple nuclear geometry either, as the 
average number of successive N+N collisions per partic- 



ipant, calculated from the Glauber model, continues to 
increase with N p in the range of the saturation. 

Presumably, the kaon yield in very peripheral heavy 
ion collisions should behave just like in N+N. As shown 
in Fig. 0, extrapolations following the trends of the K + 
and K~ yields in both Si+A and Au+Au collisions are 
consistent with the N+N values. 

Another way to characterize the dependence of kaon 
production with centrality is shown in Fig. ^ Instead 
of the average yield per participant displayed in Fig. [t], 
this figure shows the differential rate of kaon production 
with increasing N p , i.e. the slopes extracted from the 
data in Fig. [|. Values were extracted separately for each 
of the three systems. For the Si data, slopes were found 
using pairs of adjacent points, while for Au+Au data, 
slopes were obtained from straight-line fits to three data 
points in a row in order to reduce the statistical uncer- 
tainties. The lowest N p point in this figure was found 
by combining the lowest data point in Fig. with an as- 
sumed point of zero yield at zero N p . These results are 
plotted in Fig. || at the corresponding average N p val- 
ues. The shaded area shown on each plot indicates the 
range of the slope extracted from a straight line fit to 
the Si+Al data, the Si+Au data, or the Si data includ- 
ing both targets shown in Fig. ||. The kaon yields per 
participant in N+N interactions are plotted at N p = 2 
as filled points, the same as in Fig. f^. For the Si-induced 
interactions, the differential rate of kaon production rises 
rapidly and then stays roughly constant at a value consid- 
erably larger than expected from N+N. For Au+Au, the 
slope increases more slowly but continuously, and even- 
tually reaches the value found for central Si+A. In other 
words, for the same increase in N p , the increase in kaon 
yield in central Au+Au is the same as in central Si+A, 
in spite of the lower energy of the Au beam than the Si 
beam. 

Care should be taken in interpreting these results as 
demonstrating that each additional participant produces 
significantly more kaons. Increasing N p corresponds to 
decreasing impact parameter. In a simple geometric 
model, this implies that all participants have an increas- 
ing probability of multiple collisions. Also, other effects 
such as rescattering of baryons and mesons or the in- 
fluence of the nuclear medium will be stronger in more 
central collisions. Fig. ^ demonstrates how kaon produc- 
tion averaged over all participants evolves as a function 
of centrality while Fig. || shows the evolution of the dif- 
ferential rate of production. 

To investigate enhancement of strangeness production 
in heavy ion collisions relative to N+N, the heavy ion 
data in Fig. are normalized by the N+N values and 
re-plotted in Fig. ||. The filled circle at (2,1) in each 
plot represent the N+N interactions. The K + yields 
in central Si+A collisions are 3 times those in N+N, a 
result that is consistent with the previously observed 
strangeness enhancement in Si+A collisions [Q. The 
normalized K + and K~ yields per participant in cen- 
tral Au+Au collisions reach higher values than in cen- 
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tral Si+A collisions. The K + and K~ yields are en- 
hanced by factors of 4.0 ± 0.3(stat.) ± 0.5(syst.) and 
2.7 ± 0.2(stat.) ± 0.3(syst.), respectively, in the most cen- 
tral bin of the Au+Au data relative to N+N. The en- 
hancement factors are only 3.0±0.2(stat.)±0.4(syst.) and 
2.2 ± 0.1(stat.) ± 0.3(syst.) in central Si+Au collision for 
K + and K~ , respectively. The systematic errors include 
those in the yields and in N v . Errors in the estimated 
kaon yields from N+N, common to all three systems, 
are not included in the systematic errors above. The 
strangeness enhancement, therefore, is more pronounced 
in central Au+Au at 11. 1A GeV/c than in central Si+Au 
at 14. 6A GeV/c. Note that for each of the three systems, 
there is more enhancement in K + than in K~ . 
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VI. CONCLUSIONS 

In summary, E-859 and E-866 have measured kaon 
production as a function of centrality in Si+Al and 
Si+Au collisions at 14. 6A GeV/c in rapidity range 0.5 < 
y < 2.1, and in Au+Au collisions at 11. 1A GeV/c in ra- 
pidity range 0.5 < y < 1.3. In all three systems with all 
centralities, the invariant transverse mass (mi) spectra 
can be described by a single exponential in m± for both 
K + and K~; the rapidity distributions can be character- 
ized by Gaussians, with the K + distributions consistently 
broader than the K~'s. 

The kaon yields in Si+A and Au+Au collisions were 
found to have different dependencies on the total num- 
ber of participants. The kaon yields from Si+A collisions 
increase approximately linearly with the total number of 
participants, whereas the kaon yields from Au+Au colli- 
sions increase faster than linearly. For the same number 
of participants (~ 50), the kaon yields in Au+Au colli- 
sions are lower than in Si+A collisions after the beam 
energy difference in kaon production rate in N+N inter- 
actions was accounted for. On the other hand, the K + 
and K~ yields in each of the three systems have quan- 
titatively the same dependence on the total number of 
participants, differing only in magnitude by a constant 
factor, in spite of the different dN/dy distributions. 

Strangeness enhancement is observed as an increase 
in the slope of the kaon yield with the total number of 
participants as well as the yield per participant. The 
enhancement in heavy-ion collisions relative to N+N in- 
teractions appears to turn on quickly in Si+A colli- 
sions and saturate at a moderate collision centrality in 
Si+Au. The enhancement increases slowly with central- 
ity in Au+Au collisions, and in central Au+Au collisions, 
reaches and exceeds the value found in central Si+A col- 
lisions. The enhancement factor of K + production are 
3.0±0.2(stat.)±0.4(syst.) and 4.0±0.3(stat.)±0.5(syst.), 
respectively, for the most central 7% Si+Au collisions and 
the most central 4% Au+Au collisions relative to N+N 
at the correponding beam energy. 



* Present Address: Lawrence Berkeley National Labora- 
tory, Berkeley, CA 94720 

f Present Address: CERN, CH-1211, Geneve 23, Switzer- 
land 

* Present Address: IKF, August-Euler-Str. 6, D-60486 
Frankfurt, Germany 

§ Present Address: KEK National Laboratory for High En- 
ergy Physics, 1-1 Oho, Tsukuba, Ibaraki 305, Japan 

** Present Address: Oak Ridge National Laboratory, Oak 
Ridge, Tennessee 37831 

Present Address: Renaissance Technologies Corp., Stony 
Brook, New York 11790 

§§ Present Address: The Institute of Physical and Chemical 
Research (RIKEN), Saitama 351-01, Japan 

ft Present Address: Center for Nuclear Study, School of 
Science, University of Tokyo, Tanashi, Tokyo 188, Japan 

*** Present Address: Yonsei University, Seoul 120-749, Ko- 
rea 

^ Present Address: Lawrence Livermore National Labora- 
tory, Livermore, CA 94550 

Present Address: Space Science Laboratory, University 
of California, Berkeley, CA 94720 
§§§ p resen t Address: Brookhaven National Laboratory, Up- 
ton, New York 11973 



[1] L. Ahle et al. (E802 Collaboration), Phys. Rev. C 57, 
R466 (1998). 

[2] P. Koch, J. Rafelski, and W. Greiner, Phys. Lett. B 123, 
151 (1983). 

[3] P. Koch, B. Miiller, and J. Rafelski, Phys. Rep. 142, 167 
(1986). 

[4] J. Rafelski, Phys. Rep. 88, 331 (1982). 
[5] J. Rafelski and B. Miiller, Phys. Rev. Lett. 48, 1066 
(1982). 

[6] U. Heinz et al, Phys. Rev. Lett. 58, 2292 (1987). 



7 



[7] T. Abbott et al. (E802 Collaboration), Phys. Rev. Lett. 

64, 847 (1990). 
[8] T. Abbott et al. (E802 Collaboration), Phys. Lett. B 291, 

341 (1992). 

[9] R. Matiello et al, Phys. Rev. Lett. 63, 1459 (1989). 
[10] H. Sorgo et al, Phys. Lett. B 271, 37 (1991). 
[11] Y. Pang, T. J. Schlagel, and S. H. Kahana, Phys. Rev. 

Lett. 68, 2743 (1992). 
[12] B.-A. Li and C. M. Ko, Phys. Rev. C 52, 2037 (1995). 
[13] T. Abbott et al. (E802 Collaboration), Nucl. Instr. & 

Meth. A290, 41 (1990). 
[14] T. Abbott et al. (E802 Collaboration), Phys. Rev. C 50, 

1024 (1994). 

[15] L. Ahle et al. (E802 Collaboration), Phys. Rev. C 58, 
3523 (1998). 

[16] The E859 Level II trigger system, W.A. Zajc et al. (E802 
Collaboration), AIP Conf. Proc. No. 243, 1991. 

[17] F. Videbaek et al. (E802 Collaboration), Nucl. Phys. 
A590, 249c (1995). 

[18] Y. Akiba et al. (E802 Collaboration), Nucl. Phys. A610, 
139c (1996). 

[19] F. Wang, PhD thesis, Columbia University, May 1996. 
[20] B.A. Cole et al. (E802 Collaboration), Nucl. Phys. A590, 
179c (1995). 

[21] D.P. Morrison, PhD thesis, Massachusetts Institute of 

Technology, May 1994. 
[22] R. Debbe et al, Nucl. Instr. & Meth. A403, 256 (1998). 
[23] L. Ahle et al. (E802 Collaboration), Phys. Lett. B 332, 

258 (1994). 

[24] The efficiency corrections do not account for a multiplic- 
ity/centrality dependence to the tracking efficiency at a 
fixed spectrometer angle setting. Estimates of the effect 
of such a dependence indicate that it would reduce the 
most central rapidity points by at most 5% and result in 
a negligible change in the total yields. 

[25] E. Albini et al, Nucl. Phys. B84, 269 (1975). 

[26] T. Abbott et al. (E802 Collaboration), Phys. Lett. B 197, 
285 (1987). 

[27] T. Abbott et al. (E802 Collaboration), Phys. Rev. C 45, 
2933 (1992). 

[28] The average energy loss of the Au beam in the target was 
corrected for. The beam energy before the correction was 
2009 GeV. The energy loss of the Si beam in the targets 
was not corrected for. 

[29] T. Abbott et al. (E802 Collaboration), Phys. Rev. C 44, 
1611 (1991). 

[30] L. Ahle et al. (E802 Collaboration), Phys. Rev. C 59, 
2173 (1999). 

[31] B. Moskowitz (E802 Collaboration), Heavy Ion Physics 
at the AGS MIT 1993, eds. G.S.F. Stephans, S.G. Stead- 
man, W.L. Kehoe, MITLNS-2158, p. 21. 

[32] H. Fesefeldt et al, Nucl. Phys. B147, 317 (1979). 

[33] M. Gazdzicki and O. Hansen, Nucl. Phys. A528, 754 
(1991). 

[34] V. Blobel et al, Nucl. Phys. B69, 454 (1974). 

[35] There is a difference of 2% in the estimated kaon yields 
per N+N collision between Si+Al and Si+Au due to dif- 
ferent compositions of protons and neutrons in the two 
targets. The difference is well within the errors of these 
estimates and is neglected in the present study. 



8 



TABLE I. Gaussian fit results (total yield N K + and N K - , and Gaussian width a K + and a K -) of the K + and K~ dN/dy 
distributions for the Si+A data with centrality bins in calibrated target multiplicity and for the Au+Au data with centrality 
bins in calibrated zero-degree energy. The corresponding cross-sections a in mb, average zero-degree energies E in GeV, and 
average numbers of participants iV p are also listed. Quoted errors are statistical only. 



System 


Bin 


a 


E 


At 


N K + 


N K - 


cr K + 






0-210 


209 


155.4 


363 


23.7±1.6 


3.76±0.28 


0.96±0.06 


0.71±0.04 




210-309 


209 


259.5 


343 


21.2±1.4 


3.30±0.32 


0.97±0.06 


0.87±0.07 




309-488 


366 


398.3 


316 


18.2±1.2 


2.80±0.24 


1.00±0.06 


0.86±0.06 


Au+Au 


488-744 


523 


616.4 


273 


14.2±0.9 


2.15±0.18 


0.99±0.05 


0.80±0.05 




744-1099 


785 


926.4 


211 


9.34±0.56 


1.51±0.10 


0.96±0.05 


0.69±0.03 




1099-1615 


1569 


1384. 


121 


4.17±0.27 


0.64±0.05 


0.91±0.05 


0.68±0.03 




1615-2300 


1569 


1754. 


48 


l.lo±U.17 


n on i n no 

U.20±0.03 


n no _i_n i o 

U.98±0.12 


U.67±0.07 




>115 


259 


21.8 


83.9 


6.13±0.29 


1.19±0.04 


0.67±0.03 


0.56±0.02 


Si+Au 


85-115 


481 


51.7 


70.0 


5.11±0.21 


1.02±0.06 


0.75±0.03 


0.67±0.04 




50-85 


741 


149.6 


42.7 


3.17±0.22 


0.61±0.03 


0.89±0.06 


0.64±0.03 




0-50 


2222 


278.8 


17.8 


0.79±0.13 


0.12±0.01 


1.16±0.19 


0.69±0.05 




>48.5 


97 


118.3 


38.7 


2.56±0.04 


0.61±0.03 


0.82±0.01 


0.71±0.04 


Si+Al 


35.5-48.5 


180 


161.8 


32.4 


1.94±0.04 


0.44±0.04 


0.90±0.02 


0.71±0.04 




22.5-35.5 


277 


229.8 


22.4 


1.29±0.03 


0.29±0.02 


0.90±0.02 


0.73±0.04 




0-22.5 


831 


318.2 


9.5 


0.42±0.01 


0.080±0.005 


1.04±0.03 


0.78±0.05 
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FIG. 1. Calibrated and target-out subtracted zero degree energy spectrum of Au+Au interactions from E866. Centrality 
bins used in the analysis are labeled 1 through 7. 
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FIG. 2. K + (left) and K ~ (right) transverse mass distributions in Si+Au for the most central multiplicity bin. The topmost 
is absolutely normalized; others are successively divided by 10 for clarity. Data points under the curved line are obtained from 
particles identified by both the time-of-flight and Cerenkov counters. Errors shown are statistical only. Systematic errors were 
estimated to be 10%. 
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FIG. 3. K + (left) and A' - (right) transverse mass distributions in Au+Au for the most central zero-degree energy bin. The 
topmost is absolutely normalized; others are successively divided by 10 for clarity. Errors shown are statistical only. Systematic 
errors were estimated to be 10%. 
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FIG. 4. Rapidity distributions of K + (top) and if - (bottom) for Si+Au (left) and Au+Au (right) collisions. Filled-in points 
for Au+Au are data measured; open ones are data reflected about j/jvjv. Errors shown are quadratic addition of statistical 
errors and 5% systematic errors for Au+Au, and only statistical errors for Si+Au. There are 10% systematic errors in the 
absolute normalization for both Si+A and Au+Au data. 
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FIG. 5. Total K + (top) and K (bottom) yields as function of the total number of participants, N p . Errors shown are 
statistical only. Systematic error on N p is 10% for Si+A data and 8% for Au+Au data. Systematic errors on the yields are 
10% for both Si+A and Au+Au data. Note that, for clarity, the ordinate scales start from negative values. 
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FIG. 6. Ratio of K + to K total yield as function of the total number of participants, N p . Errors shown are statistical only. 
Systematic error on N p is 10% for Si+A data and 8% for Au+Au data. 
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FIG. 7. Total K + (top) and K (bottom) yields per participant as function of the total number of participants, N p . Errors 
shown are statistical only. Systematic error on N p is 10% for Si+A data and 8% for Au+Au data. Systematic errors on the 
yields are 10% for both Si+A and Au+Au data. Filled points are those estimated for N+N interactions at y/s — 4.74 GeV 
(filled circles) and 5.39 GeV (filled squares). 
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FIG. 8. Differential production rate for K + (top) and K~ (bottom) as function of the total number of participants, 7V P . 
Errors shown are statistical only. Filled points are the kaon yields per participant in N+N interactions at y/s = 4.74 GeV (filled 
circles) and 5.39 GeV (filled squares). See text for explanation of the shaded areas and details on the method of extracting the 
differential production rate. 
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FIG. 9. Total if 4 " (top) and K (bottom) yields per participant, normalized by those in N+N interactions, as function of 
the total number of participants, N p . Errors shown are statistical only. See text for discussion of the systematic errors. 
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